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ditions is calculated exactly as a series in the absence of an external magnetic field. 
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'Th 'I I. INTRODUCTION 

o ; 

The free energy of conformally invariant two dimensional systems have been studied for some time [l| 0] • In the 
case where the system considered is an infinite strip, it is possible to find an exact expression for the free enery per 
unit length. This has been done in the isotropic case, with equal horizontal and vertical coupling constants [i] . This 
was done by considering a cylindrical Ising lattice of height Ai and circumference 2Af with so called Brascamp-Kunz 
boundary conditions [5[ . In this paper, the corresponding calculation in the case of different vertical and horizontal 
coupling constants will be done. 

Blote, Cardy and Nightingale [l| wrote the limit of the free energy at criticality per unit length limw-too F/2AT as 



where / is the bulk free energy per unit area and / x /2 is the surface free energy. In the isotropic case A is given by 

A = --(c-24/i) (2) 

where c — 24ft, is called the effective central charge. The effective central charge has been discussed elsewhere, for 
\Q ' instance by Izmailian, Priezzhev, Ruelle and Hu 6]. For the Ising lattice, the central charge is c = 1/2, and the 
I/"") , allowed values of the conformal weight h are 0, 1/2 and 1/16. It can be seen from equation (17) of ref. 0] that 
the partition function is a multiple of the Virasoro character Xi/w — y/fa/^V (where the suppressed argument is 
iN /(M + 1)). Thus h must be 1/16. It will be shown later in this paper that in the isotropic case 

A=£. (3) 

O ■ Since c = 1/2, © and © confirm that h = 1/16. 

J> \ Equation ([3]) is only valid if the classical system is rotationally invariant at large distances [H 0- I n terms of 
an Ising lattice, this means that the lattice has to be isotropic. If this is not the case, then ^ must be modified 
by dividing the right hand side by v, the "speed of light". The speed of light should be obtained from a dispersion 
^ ■ relation uj ~ vu, where u> is the frequency and u is the momentum. 

The two dimensional cylindrical Ising lattice introduced by Brascamp and Kunz 0] is the lattice A = {(m, n) | 1 < 
77i < Ai, 1 < n < 2Af, (m, 2Af + 1) = (m, 1)} with the following boundary conditions: 

(i) The lattice interacts with a row of fixed, positive spins above it, 

(ii) The lattice interacts with a row of fixed, alternating spins below it. 

(See fig. [l}. The interaction between neighboring spins on A is E\ in the horizontal direction and E% in the vertical 
direction. In what follows, the dimensionless parameters Ki — j3Ei (1 = 1, 2) will be used instead of E\ and E^. 
Apart from from interactions between nearest neighbors, there is an external magnetic field. Thus the Hamiltonian 
is defined as 

M 1M M 2M M 2M 

£ A {<7, Kx,K 2 , H) = -Ei ^ X! ~ E 2 X! X] a J,k<yj+i,k - y^^~] H(j, k)<jj, k (4) 

j=l fc=l j=0 k=l j=l k=l 
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FIG. 1: The lattice A (marked by "x"; A4 — 3, J\f = 2) and the dual lattice A* (at intersections of lines). 



where ery^ = <Jj t k+2N = ilj a o,k — 1 arid oa-i+i./c = (— 1) The partition function 



Z A (Jfi ,K 2 ,H) = ex P ( CT > > ^ , # ) 

cre{-l,l} A 



(5) 



has been calculated for the constant external magnetic fields H = and = in/ 2 for which the problem is exactly 
solvable. Brascamp and Kunz calculated Z\(K,K, 0). Z A {K\, K 2 , 0) and Z A {K\, K 2 , in/2) were calculated in ref. 
0. Z A (#i,# 2 ,0) is given by 

Z A (Ki,K 2 ,0) = 2 2MAf J | J| {cosh 2^1 cosh 2^2 - sinh 2 JTi cos 6j - sinh2if 2 cos^ fe } (6) 
i=i k=i 

where 



(2j - I)vr/2AA, tp k = kn/{M + 1). 



(7) 



II. THE FREE ENERGY 



The partition function ((6J can be written as 



Z A = 2 M "e 2M ^ K ^ ]J H F(j, k), 



TV M 

nn 

3 = 1 k=l 



(8) 



where 



/ / , 9 /sinh 2K 2 x 1/2 

F ^ : = 4 2sinh U*2j£ 



i, #2) + sin 2 0,-/2 + sin 2 ^/2 

smh2Ai 



(9) 



and the mass /J,(Ki, K 2 ) is defined as 



sinh 2 fi(Ki,K 2 ) 



4(8^2^ sinh 2K 2 y/ 



jj5 (cosh 2Ki cosh 2K 2 - sinh2A'i - sinh 2K 2 ); 



in particular 



(10) 



H{K,K) = - log sinh 2K. 



(11) 
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, n , . , n ( sinh^A^X , \ sinh2AT 2 9 , n 

sinh 2 u {Ki, K 2] u) := 2 sinh 2 — — 2 - ^K l3 K 2 ) + -— # sin 2 u. (12) 



Define u(Ki, \ u) by the equation 

Vsinh2ifiy' ^ v ' 1 sinh2Xi 
so that 

F(j,k) = 4(sinh 2 w(^i,X 2 ;^/2)+sin 2 ^/2). (13) 

Then, using the identity Q 

2JV-1 

JJ 4 (sinh 2 w + sin 2 6^/2) = 4 cosh 2 2A/w (14) 
i=o 

one obtains 

2AT-1 

JJ Ftf + l,0)F(j + l,X + l) = [4cosh2AAw(^ 1 ,if 2 ;0)cosh2^(is: 1 ,^ 2 ;7r/2)] 2 . (15) 
3=0 

Therefore, using the same argument as in ref. (4j , one finds 

22M/V e 4MA0i(Ki,Ki) 

Zl = t-, ; r^X (1/2, 0, n(K 1 ,K 2 ) ) (16) 

A 4cosh2AM^i,^2;0)cosh2A/^(A:i,A: 2 ;7r/2) 

where A is a lattice of size (2(vW + 1), 2Af) and the partition function of A with twisted boundary conditions is given 
by the equation 

2AT-1 2M+1 



( z-„ M - n n <-«^w^ + ^^^). <"> 

(Z A (l/2,0, Ai (A' 1 ,^ 2 ))) 2 is given by 



(Z A (l/2,0 J /i(Jfi,Jf 2 ))) a = II II 4(2 sinh 2 /i(A-!, g 2 ) + sin 2 g J+1 /2+ ^ h2j ^ 2 sin 2 y fc /2 

2JV-1 2M+1 

= n n ni+^k) 

j=0 fc=0 



sinh2ATi 

2^-1 2M + 1 f , . x 1/2 



TT TT ( ( — I (cosh 2ATi cosh 2AT 2 - sinh 2K X - sinh 2AT 2 ) 

x -*- A - 1 - V V sinh 2 A" 2 / 

j=0 k=Q v v 7 

+2 sinh 2 ATi sin 2 9 j+1 /2 + 2 sinh 2 AT 2 sin 2 p k /2 ] . (18) 



Define Q(Ki, A" 2 ; u) by the lattice dispersion relation 

sinh 2 ^(A' 1 ,X 2 ;u) - ^ (2 sinh 2 ii(K x ,K 2 )+ sin 2 u) 
sinh 2A 2 



(sinh 2Xl)1/2 (cosh 2A^ cosh 2A- 2 - sinh 2A"!- sinh 2A- 2 ) + ^^i sin 2 (19) 



2 (sinh 2 AT 2 ) 3/2 sinh2A" ; 
At criticality it reads 



One would thus expect that 

/sinh2ifi x 1/2 



v 



\sinh 2K 2 

It will be shown later that this is in fact the case. It follows from (|18p and (fT9|) that 

4Af(M+l) 2AT-1 2M+1 



^y. 2g 2 ) II II 4(sinh 2 a;(ir 1 , j pr 2 ;^ +1 /2) + S in 

1 ' j=o k=o 

Using the identity 

2M+1 

Yl 4 (sinh 2 lu + sin 2 tp k /2) = 4sinh 2 2(A1 + l)w 



2 



fc=0 

one obtains from ([22]) 

2AT(M+1) 2AT-1 



^ h2 ^ 2 ) II 2smh2(M + l)a}(iri,ir 2 ;^/2) 

Mill A j. / 



III. ASYMPTOTIC EXPANSION OF THE FREE ENERGY 



Let (Ki,K 2 ) = (K* ,K) be the curve on which ^(Ki,K 2 ) = 0, or equivalently 

sinh 2.^1 sinh 2^2 = 1. 

Let the Taylor expansion of oj(K* , K; u) be 

A2n 



,2n+l 



-** (2n)\ 

n=0 v ' 



In particular 



Then 



and similarly 



/ sinh 2K*\ 1/2 . 

A = = sinh 2K* 

V sinh2K J 



lv(K* , K;u) = sinh 1 Aq sin u = sinh 1 sinh 2K* sin u 



u(K* , K; u) = sinh 1 A 1 sin u = sinh 1 sinh 2K sin u. 
Thus, on the critical curve, the free energy as obtained from (116[) is 

F = -logZ A = -TWA/log 2 - 2MMn(K*,K*) + log2 

1 /2 

+ 1 log cosh 27V" sinh' 1 f j - \ log Z A (l/2, 0, fi(K*, K)) 

= -MAflog2-2MAfp,(K*,K*) + log2 

+ X - log cosh AMK - l - log Z A (l/2, 0, n(K*, K)). 

Clearly 

logZ A (l/2,0, M (X*,if)) = 2M(M + 1) log ( * + 2(A4 + 1) ]T a)(^*,iC,^/2) 

2AT-1 

+ ^ log(l-exp-4(7W+l)a}(Jf*,ir,^/2)). 
j=o 

The two sums in (f5T|) can be calculated exactly up to an exponentially small correction 0{e~^). 
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A. Calculation of ((31]) 

The first sum in (|3ip can be written as a power series using the Euler-Maclaurin summation formula: 



2A ^ ±1 ^ — A 2 „ B 2n+2 (l/2) 



2(M + 1) V D(K*,K,9 j /2) = - / w (if*, u) cfct - 2tt£ V 

j=0 ^ ^ n=0 



(2n)! 2n + 2 



(32) 



where S* = 4Af(M + 1), £ = (M + 1)/N and B fe (l/2) is the Bernoulli function B k (x) evaluated at x = 1/2. It is 
defined as 

B p (a;) := V k- p e 27Tlkx . (33) 

v ; fcez\{o} 

It remains to calculate the second sum of (f3Tj) . 
The second sum in (1311) can be written as 



2AT-1 oo Af-1 



E log (l - e -i(M+l)^(K%K,6 J /2)^ = -2 E — E e - 2m2{M+1)ij{K ' '> K ' 6 i/ 2 ') . (34) 



m 

j—0 m— 1 j— 



Let P(p) = {7r = (qi, g„, ri, r v ) q 3 ,rj e N, 1 < ^ < p, ^ ^ g fe if j ^ k, 1j r j = P}- Tne exponential 

the right hand side of (fM)) can be written as 



-2m2(M + l)Cb{K* ,K,6 j/2) 



= exp (-2nm\ Z(j + 1/2) - 2nm^ E (j + l/2) 2 ?+^ 

= (l - a^f) (J+ ( 1 2 / J P+1 A 2P ) c-^tf+Va) (35) 



where 



A 2P = (2p)! E I II i (t^T ) ] ("2™f (? + l/2)r + -^W" 1 (36) 



^eP(p) \i=i ' v ^ y(y 
Together, {34]), (J35J) and (J3BJ imply that 

2JV-1 oo /jV-1 

E lo S (l - e-^ 1 ^.^^) = -2 E ^ E 

3=0 m=l \ 3 =0 



e -27rmAoC(j+l/2) 



00 / 2f\P i / oo A/"— 1 \ 

+ 4 < E ( V) ^ EE^ + 1/2) 2 P+ l e -2^A f(3 + I/2) 1 (37) 

As in ref. if, for large J\T, the finite sum Y^Lq 1 ^ s replaced by the infinite sum YlfLo m (ETJ) , then equality still 
holds up to an exponentially small correction Oie^^). Thus 

2Af— 1 oo / oo 

-27rmA f(j+l/2) 



E lQ g (l - e-^W*^) = 2 E ^ E ' 

j=0 m = l Vj=0 



4 <E(^) P (^W( EE^^re^^l+Ote^). (38) 



OO OO 
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Combining (f3"2")l . and (f5T|) . one obtains 

log %(l/2, 0,//(ir, if)) - 2TV(.M + 1) log ( jgj^ 1 + Z / u(K*,K,u)du 



Vsinh 2/sT* 

n=0 x 7 v y m=l j=0 

00 / 2 £ \ P i oo oc 

+ ke t ^jj;(,-+i/2r e -^v«) +0(e ^ ) ( 39) 

This expression may be further simplified in terms of elliptic theta functions. However, it is simplest to consider the 
limit TV — > oo. 



B. The free energy in the limit TV —* oo 

In the limit TV — ^ oo while TW is fixed, an exact result can be obtained. 
Combining (|3"0|) and (f3"9")l one obtains 

= -~Mlog2-(M + 2)n{K,K)+K 



_ l im ___±1 r Q(K*,K,u)du+- hm 1 V 1 V e ^mA eO + i/2) 

Af-oo f L 1 ' ' ; 2 AA^oo TV ^ TO ^ 

J ° m=l j=0 

i 00 / 2 \ P 1 oooo 

- JfiL i<E^ + l/2)«'«— «««« . (40) 



p=l v / \ V) m= i j = o 

The limit of the double sum can be calculated to be 

1 ^ oo * oo 

2 A^™oo TV ^ m ^ 24(X + 1)A [ ' 

m—1 3=0 

Further, by ref. @ 

Y: £(j + i/^e-^W+i/a) = _^(B 2p+2 (l/2) - K^i*,*)) (42) 

m=l 3=0 ^ ' 

where B p (x) is the pth Bernoulli function and 

| -2iri(na+mf3) 

K a '"(T) •= V (43) 

p K ' (-2m)P ^ (n + rm)P K ' 

(m,n)^(0,0) 

is Kronecker's double series. 

It can be shown 01 that K 2 p 2,0 (i£) can be expressed in terms of the elliptic theta functions $ 2 , $3 and $4. It can 
therefore be shown [J] that for small £ 

k^aoO = B 2p+2 (A er 2p - 2 + ocr*- 1 ) (44) 

where B„ := B„(l) is the nth Bernoulli number. Since 

A 2p = X 2p + 0(0 (45) 
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for small £, it follws that 

f i m + 1 r 

K 



24{M + 1)A 



2p+l o 

D9 



+ V - (46) 



In particular, one sees that the expected value of v given in (|2"Tj) is correct. 



C. The free energy for large J\f 



We now consider the case of large A/", where the partition function is given by (|39|) . According to ref. [9(, 



OG -. OC 



_ 2 V i f e -^A D fl J+ i/2) = 2 V log (1 - e - 2 ^^^+ 1 /2) ) = log ^(»AoC) + 7 rA gB,(l/2) (47) 

where r)(r) := ($2(t)M t )M t )/ 2 ) 1/3 = e l7TT / 12 ]J^ =1 {l - e i27rrn ). From (39]), (J42J) and 07]) one thus obtains 
logZ A (l/2,Q,KK*,K)) = 2N{M + 1) log ^j+^£u(K*,K, u) du 

£z£ V s J i 2n ) 1 2n + 2 ?7(«AoC) 

+ <E (V) (2^)! A2p ^TT (B2p+2(1/2) " K ^ A °^ + °( e_A ( 48 ) 
It now follows from and (gSJ) that 



F = -.MAnog2-2A^A/>t(A*,A*) + ilog2 

/ sinh 2 A \ S r 



<D (A*, A, u) <At 



^ , ^A" A 2n B 2n+2 (i/2) i , , fl 4 (a fl iw Rn/9 s 

n=0 x x ' 

2<E (V) (2rf! A2p ^TT (B2p+2(1/2) ~ K ^+° (iA ° 0) + ° (e ^ } (49) 
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